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Abstract 

Different conclusions about quark nature of /o(980) based on the analysis of data on 
the reaction c/>(1020) — > 7/0 (980) make it necessary to perform detailed comparison of 
formulae used by various groups for the calculation of partial widths of radiative decays. 
We carry out a comparative analysis of methods of calculation of radiative decays 
4> -» 7/0 and -> 777 performed by F.E. Close et al, Phys. Rev. D 65, 092003 (201 
and A.V. Anisovich et al, Eur. Phys. J. A 12, 103 (2001); Yad. Fiz. 65, 523 (2002). 
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We carry out a comparative analysis of methods of calculation of radiative decays like 
<f> -» 7/0 and (f> -> 777 performed by F.E. Close et al, Phys. Rev. D 65, 092003 (2002), 

Q-t, 

^ ■ 1 Introduction 

> : 

At present there is a common understanding that radiative decays are a powerful tool for the 
study of the quark structure of mesons, and the calculation of corresponding amplitudes is 
a subject of the increasing interest. Several groups are involved in such calculations, with 
different results and different conclusions about the quark structure of mesons based on these 
calculations. One of the main differences concerns the quark structure of / (980) and a (980). 
According to [UI2], /o(980) and a (980) are dominantly qqqq states, and the authors claimed the 
hypothesis about a dominant qq component to be contradicting the data on radiative decays. 
At the same time, in the papers of our group [13 lU E], an opposite conclusion is made: the 
widths of radiative decays agree with the hypothesis about the qq origin of / (980) and a (980), 
and the ss component in / (980) should be about 60%. The sharp discrepancy of conclusions 
demands to clarify the points, which resulted in such a disagreement. 

The recent paper jHj, where the calculation of a variety of reactions of one-photon radiative 
decays was carried out, provided us with an opportunity to compare the methods of calculation. 
Still, in jS], the calculation was performed in non-relativistic approximation, while in jUllHE] we 
adopted a relativistic approach for the gg-systems. Let us point out that in j3J |3] , the formulae 
for the amplitudes of radiative decays were written in a form allowing a simple transformation to 
non-relativistic approximation for quark composite systems. The result of comparison appeared 
to be startling: in the transition amplitudes depend explicitly on the momentum of the 
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decay particles, while in |U E] there is no such dependence. This circumstamce enforces us to 
perform a comparative analysis of the two methods — this is the subject of the present paper. 

We consider form factors of radiative decays involving vector (V), scalar (S) and pseu- 
doscalar (P) mesons, such as V — > 7P and V — > ^S, written in terms of Feynman integrals; 
then we transform them to non-relativistic approxomation. In the framework of such a proce- 
dure, one can see directly that in the amplitude of radiative decays V — > 7P and V —* jS the 
momentum depependence, like that obtained in |6., is absent. The momentum dependence of 
the matrix element implemented in [U] causes the violation of Lorentz-invariance: the require- 
ment of Lorentz-invariance, due to zero photon mass, reveals itself also in a non-relativistic 
description of composite systems. 

Let us clarify this statement in more detail. The Lorentz-invariant structure of the ampli- 
tudes of the transitions V — > 7P and V — > '-/S have the form: 

Ay^p = -ie^^e^e^q^p^Fv^piq 2 ) , 

Here e^ auiV2 is the antisymmetrical tensor, and are polarization four-vectors for the 
vector meson and photon, and p and q are their four-momenta; g m is the metric tensor. Invari- 
ant form factor F v ^ lP (q 2 ) (or Py^ 7 s(g 2 )) is a function of three independent variables: m v , q 2 
and m 2 p (or m§). For real decays q 2 = 0, and this circumstance simplifies calculations of the 
form factors Fy^ jP (0) and Py^ 7 ,s(0) as convolutions of quark wave functions. Form factors 
represented in terms of light-cone variables k±_ and x read as follows [3J |5] : 

Fy^ P (0) ~ ip v (k ± ,x) <g> ip P (k ± ,x) , 

F v ^ s (0) ~ ifr(k x ,x) <g> if>s(k±,x) , (2) 

and in these convolutions the photon momentum q is not present explicitly. In the relativistic 
form factor representation, such a property is obvious because of the Lorentz invariance: form 
factors with q 2 = allow one to choose any set of variables, go = |q| = included. 

After coming to non-relativistic approach, the form factors obtained in |3J |5] turn into 
the convolutions of quark wave functions: 

^i(|k|) ® ^(|k|) , (3) 

where k is relative momentum of the gg-system. Another representation of form factors is 
declared in [S|: form factors for the transitions V — > 7P and V —> 'jS are written as convolutions 



^i(|k|) ® i) 2 



k+^q 



(4) 



(to avoid misunderstanding, let us note that in [Hj the photon momentum was denoted as p, 
see Eq. (2) of this paper, but here we use q for it). 

The formulae and (jlj reflect the discrepancy between papers |3J EJ |3] and jB] , so they 
are the subject of the present discussion. 
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By discussing non-relativistic formulae for the transitions V — > 7P and V — > jS, it would 
be more transparent not to use the trick with the choice of variables but to work directly with 
physical momenta qo = |q| 7^ 0. Such a calculation is done in the next Section by studying 
form factors Fy^ 7 p(0) and Fv->ys(fy i n the non-relativistic limit. We show that, by making 
a correct transformation, when in the integrand the terms of different order are not kept, the 
radiative decay form factors are defined by the convolution 

Another difference between papers [6] and [31 EJ [5] consists in a way of the application of 
the threshold theorem [7] (i*V^ 7 s(0) ~ uj at small mass difference of the vector and scalar 
particles, uj = m v — m s ) to the decay 0(1020) — > 7/0 (980). In [3J HJ [5] the process V — > jS 
was considered in terms of the additive quark model when the photon is emitted by one of the 
constituent quarks. For this process the equality Fv~^s(0) = is impossible if V and S are 
basic states (l 3 Siqq and l 3 Poqq), for radial wave functions of these states have no nodes and 
the convolution J5J) does not change sign. For basic states the requirement Fy^ 7 s(0) = is 
fulfilled only under the inclusion of processes beyond the additive model jH], for example, it can 
be the photon emission by the t-channel exchange current [§|. In jH] the threshold behaviour, 
Fy^ 7 s(0) ~ uj at uj — > (see Eq. (10) of paper 0), is believed to be due to the proper choice 
of the wave functions. 

We do not investigate here the threshold behaviour of form factors in detail - special articles 
are devoted to this subject [1113, where two types of reactions have been considered: with the 
production of a bound system and with production of two unbound particles with the same 
quantum numbers, for example, 0(1020) — > 7/ (980) and e + e~ — > —* 7(^)5. The amplitudes 
of these reactions are connected with each other, namely, the transition amplitude of 0(1020) — > 
7/o(980) is a residue in the poles of the amplitude of the process e + e~ — > (f> — > 7(^^)5: 



(0) 

ml - s e+e -)(m} 



A e + e -^ 9!) _ +7 ( 7r7r ) s (s e + e - ,s n7r ,0) = -. — 5 — \/° 2 \ + smo °ther contributions , (5) 



where s e + e - and s n7r are the effective masses squared of the e + e and im systems. The threshold 
theorem is valid for the reactions with stable particles, that is for A e + e -_^.^_» 7 ( 7r7r ) s (s e + e - , s nn , 0) 
and ^40(io2o)-»7(7T7r) s ( m % s ttit> 0): 

^■e+e - — >0— >7(7T7r)s ("^e+e - 5 3-ir-m 0) ~ (s e + e - ^7T7r) 5 (6) 



^V(1020H7(7T7r) s ( m ^, Stttt, 0) ~ (m^ 



m,-s 



7V7V ) •} 



for the 0(1020) can be considered as a stable particle with a good accuracy. 

The problem of direct application of the threshold theorem to the /o(980) is not self-evident 
because this resonance definitely cannot be considered as a stable particle: it is characterized 
by two poles on different sheets related to the KK threshold, at 1020 — i40 MeV and 960 — «200 
MeV, and the second pole is essential for the formation of the /o(980) signal, see 9\ for more 
detail. 
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2 Form factors i*V_> 7 s(0) an d p~V^ 7 p(0) m additive quark 
model 



In the additive quark model, radiative decay is a three-stage process: the transition V — > gg, 
photon emission by one of the quarks and the fusion of quarks into a final meson (S or P), 
see Fig.l. The processes considered here, V — > 75 and V — > 7P, are the transitions of both 
electric and magnetic types. Accordingly, the formula for the photon-quark vertex reads: 

where m is the quark mass, a a p = (7«7/3 — 7/?7a)/2, for the notations of momenta see Fig. 
1. Such a representation of the vertex being eqivalent to the expression with the use of 7^ 
simplifies the calculations related to the transform to non-relativistic limit. 

In the calculations, we work with amplitudues written as e^~ , e^A^" ,p ^ s taking into ac- 
count the requirements e^q a = and e^p^ = 0. Therefore, the calculated amplitudes obey 
the constraints q a AY^* r,p / s = and p^A^" ,p / s = 0. 



2.1 Feynman representation of the triangle diagram of Fig. 1 

Feynman integral for the diagram of Fig. 1 reads: 



1 - 1 1 \ 1 / d 4 k „ (-)Sp 7m0i + m)T a (k' 1 +m)T(-k 2 + m) 

4^ ' I -TTT^jGv- 



G , (8) 



z(2tt) 4 v (m 2 - kl - i0)(m 2 - kf - i0)(m 2 -k\- zO) 
where for the transition V — ► 7S the spin operators T and T a are defined as follows: 

r = /, T a = ±-{k la + k' la) , (9) 

and for V — > 7P they are equal to 

F = 75 , F a = cr a/37 -^- . (10) 

2m 

Transition vertices V — > gg and qq —>■ S (or gg — ► P) are denoted as GV an d G. 

A suitable tranformation procedure for getting non-relativistic expression is to introduce 
in (JBJ) the two-component spinors for quark and antiquark, ipj and Xji that is realized by 
susbstituting 

( kl + m ) ^ £ u j (ki)u j (kx) , ft + m) - J] ^'(Wft) , (11) 

5=1,2 j=l,2 



where 

u j (k) = \/k + m \ 

ko+m 



= vA ., ,„ I ^ J . ./= 1.2. (12) 



and 

(k-m) - u\-k 2 )u\-k 2 ) , (13) 

i=3,4 

where 

u^—k) = i^/k + m ■ ( 14 ) 

After substitutions pij) -(|14 j) . we have the two-dimensional trace in the numerator of the inte- 
grand (jHJ). 



2.2 The diagram of Fig.l in non-relativistic approximation 

Now we go to non-relativistic approximation in the vector-particle rest frame. Denoting the 
four-momentum of vector particle as p = (po, p_i_,p 2 ), we have in this frame: 

p = (2m-e v ,0,0) , (15) 

where Ey is the binding energy of vector particle, which is supposed to be small as compared 
to the quark mass, Ey <C m. 

Let the photon fly along the z-axis, then 

q = (q z ,0,q z ) , (16) 

and the polarization vector of photon lays in the (x, ?/)-plane. 

The four-momentum of scalar (pseudoscalar) particle is equal to 

p' = (2m -e v - q z , 0, -q z ) = (2m - e) 2 + ql 0, -^2 m - £ + |-, 0, -q)j . (17) 

Here e is the binding energy of scalar (pseudoscalar) particle, which is also small compared to 
the mass of constituent e <C m . 



2.3 The reaction V ^ 

The transition to non-relativistic approximation in the numerator of the integrand (JHJ) provides 
the following formula for the reaction V — > ^S: 

- Sp 2 [2m^ ■ (k 2 - k'Jo-] (k la + k' la ) . (18) 

The notation Sp2 stands for the trace of two-dimensional matrices. In (j!8|) . in the transition 
to non-relativism, the following terms are of the leading order. 
In the photon-quark vertex (electric interaction): 

u(*i) ha 2 m k[a — > V + (mi, + k'MV) , (19) 



in the vertex V — > qq: 

H-hh^uih) — > x + (2)2ma^(l) , (20) 

in the vertex qq — > S 1 : 

— v^ + (l / )^(k 2 -k' 1 ) X (2) . (21) 

The constituent propagators in the transition to non-relativism should be replaced as follows: 

(m 2 - k 2 - iO)- 1 — ► {-2mE + k 2 - iO)' 1 , (22) 

where E = ko — m and m 2 — &q ~ —2mE. Then the amplitude for the transition V — > 7^ is 
defined by the diagram of Fig. 1, it reads as follows: 

(V) ( 7 ) f dEd3k r -Sp 2 [2mff M ■ (kg - k / i)o-](fci Q , + k' la ) 

" Cq 7 i(2vr) 4 y (-2m£?i + k 2 - iO)(-2m£J + kf - tO)(-2mE 2 + k 2 - iO) { } 

Furthermore, we denote E 2 = -E, k 2 = k. With these notations the energy-momentum conse- 
vation laws are written as follows: 

Ei = —sy — E , ki = — k , 
E[ = -e v - E - q z , ki = -k - q . (24) 

One can integrate over E by closing the integration contour in the lower half-plane of the 
complex- valued E, that is equivalent to substitution in (J2H|) : 

( _ 2mB+k2 _ iorl _^( £ _jg. (25) 

By fixing E = k 2 /2m, we can evaluate the order of value of momenta entering (|23|). Formula 
(|17|) gives us 2m — ey — q z = 2m — e + q 2 /2m 2 , or 

g z ~ £y - e , (26) 

because q 2 z /2m is the magnitude of the next-to-leading order of value. In the integrand of (J23J), 
the momentum square is of the order of 

k 2 ~ 2me ~ 2me v , (27) 

and this means that 

Qz < |k| . (28) 

Therefore, within non-relativistic approximation, the amplitude for the transition V — > 7S 
reads: 

r H 3 k 1 
e ^ 7) / ^^ Ws(fc) 2^ ( ~ 4) Sp 2 [2m^-2(fc<r)](-2A; Q ) , (29) 



where the requirement (j2%|) is duely taken into account and the wave functions for vector and 
scalar particle are introduced: 

Mk) = T f 1 5\ > = Al ~, , • (3°) 

A(me v + k 2 ) v ' 4(me + k 2 ) v y 



Recall, polarization vector do not contain the time-like component and polarization vector 
of the photon belongs to the (x, y)-plane. Accounting for Sp 2 [c^crg] = 25^, where 5^ is the 
three-dimensional Kronecker symbol, and substituting in the integrand k^ka — > Sfj,ak 2 /3, we 
have had final expression: 

Ay^s=(e {V) e^) J -j^Mk) M*>) f *? = / ^ MQMQ^ * 3 . (31) 

\ ) •■■> 

Here we re-denoted k 2 — > /c 2 . The photon polarization belongs to the plane orthogonal to p 
and g, so one can re- write (|31|) in the form of Eq. (0) by using 

K m ipq) J " ~ 

Formula ()31|) coincides with that obtained in terms of the double spectral-integral representation 
(Eq. (32) of =4]), if in the relativistic- invariant formula (32) of [3] one makes a transform to the 
non-relativistic limit. 



2.4 The reaction V — > 7P 

In the reaction V —>■ jP, non-relativistic spin factor (numerator of the integrand of (jHJ) has the 
form: 

(-) Sp 2 [2ma M ■ ie af3y q^a 1 ■ 2m] = -iSrrfe^p q p , (32) 
where e a ^ is three-dimensional antisymmetric tensor. As a result, we have: 

^ v {k)t\) P {k) . (33) 

7T 7T 



This expression coincides also with its relativistic-invariant counterpart obtained in |2] (Eq. 
(61)), if in |4] one makes a transition to non-relativistic limit. 



2.5 Normalization of wave functions ipy(k), ips{k) and ipp{k) 

Here we use the wave function normalization similar to what was applied in [31 El E] for the 
relativistic treatment of composite systems. 

For scalar particle, the normalization reads: 



and for vector particle it is 



/ — mr— = 1 , (34) 

J IT TTm 



d ^L tfr(k)— = 1 . (35) 

7T 7T 



Normalization for pseudoscalr particle is the same as for vector one: in (J35|) one should replace 
Vv(fc) -> 4>p(k). 

Normalization condition of wave function can be re-formulated as requirement for charge 
form factor at q 2 = 0: 

F charge (0) = 1 • (36) 

It is easy to see that formulae (JMj) an d (|35p are actually the requirements for charge form 
factors. Let us explain this using scalar particle as an example. 

Form factor of scalar particle is defined by triangle diagram of Fig. 1 type. Using the same 
calculation technique, which resulted in formula (|29|1 . we obtain: 

/d^h 2 1 

^^{k)-$ V2 [2{k*).2{ka)\- . (37) 

Here, as by derivation of Eq. (|29|) . the factor 2/m arises due to the integration over E and 
wave function definition ips(k); the vertex S — > qq is equal to 2(fcer), and the factor 1/2 
appeared because of substitution k\ a + k' la — >■ (pi Q + j/ lc J/2: at a = 0, that corresponds to 
the interaction with Coulomb field, we have k\o — k' l0 ~ m and pio = p' 10 — 2m). Just the 
condition (0) = 1 gives us the formula (JHlj) . 



3 Comments on the analyticity of the amplitude V — > 7*S 

1. The threshold theorem can be reformulated as the requirement of analyticity of the field 
theory amplitude. The transition amplitude V — ► 76* written for elementary particles of the 
field theory Lagrangian has a structure as follows: 

g, a -^^-)A(s v ,s s ,0) , (38) 



sv - s s , 

where (pq) = (sy — ss)/2 and g 2 = 0. Then analyticity constraint reads 

[A(s v ,s s ,0)) sv _> ss ^0. (39) 

The constraint (|39j) is written for the particles which form the standard sets of \in) and (out\ 
states, i.e. for particles which can be considered as stable ones. The problem of applying (|3*U|) to 
resonances is questionable because the resonance amplitude is determined by the pole residue, 
while the analytisity requirement (0) includes smoother terms as well. 

2. The amplitude of additive quark model for the transition V — > determined by (|31|) 
cannot be equal to zero. Indeed, to turn the right-hand side of Eq. (31) to zero one would 
need at least that one wave function was a sign-changing one. But wave functions of the lowest 
radial-excitation states, with radial quantum number n = 1, have no nodes. Therefore, the 
transition matrix element for basic states with n = 1 does not turn to zero, this fact not 
depending on what potential for V and S systems is under consideration. 
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Although in the additive approach Ay^sify 7^ 0, the factor (my — ms) inherent in the 
El-transition can be extracted explicitly from the right-hand side of (|31|) for a wide class of 
models. 

In order to clarify this point let us consider as an example the exponetial approximation for 
the wave functions ipy(k 2 ) and ips{k 2 ): 

ipv{q) ~ exp , i>s(q) ~ exp (-6s A; 2 ) • (40) 

With exponential wave functions, the matrix element for the V — > amplitude given by the 
additive quark model diagram, Eq. (|3ip. up to a numerical factor is equal to: 

,3/4,5/4 

Ay^s ~ (e e ) ( ^ + ^ )5/2 - (41) 

In case of one-flavour quarks in the exponetial potential well one has 

m(my — ms) = by 1 , (42) 

that gives the amplitude Ay^$ i n the dipole emission representation: 

,7/4,5/4 

Ay^s ~ (W)) ^, 5 , 5/2 m(m v - m s ) . (43) 

In case under consideration the factor (my —ms), or (£5 — 5y), relates to the difference between 
the V and 5 levels and is defined by the by value only The amplitude v4y^ 7 s(0) turns to zero 
when binding energies of either ^/-meson or S'-meson tend to zero: Ey — > or 63 — > 0. 

The considered example demonstrates that the statement of Ref. that the triangle 
diagram contribution does not contain the factor (eg — By) is not generally justified. Still, the 
equivalence of the additive model and dipole representations for transition amplitudes depends 
on the type of interaction under consideration, see jU] for more detail. 

3. There is one more point that is worth mentioning in connection with the present discus- 
sion, that is a representation of the spin operator. In 0, it was emphasized that the choice of 
the spin operator (J3SJ) is not unique as sometimes was it suggested [T]. For the spin structure 
of the amplitude V —> 7S' with emission of a real photon (q 2 = 0), one can alternatively use 
the metric tensor 5^(0) defined in the space orthogonal to the four-momenta p and q: 



#ua(°) = Qua + , s 2 gpga 1 (jPuQa + QfiPa) • (44) 

[sv - s s ) 2 sv- s s 

The matter is that the difference 

^(0) - (gya - P&L. ) = 4L^(0) , ( Jo) 



where 



s v - s s 



L^(0) = SV q^qa - — ?P»qa, (46) 

(Sv-Ss) 2 2(S V -S S ) 
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is a nilpotent operator: 

W n ) W°) = . (47) 

Adding a nilpotent operator to a spin operator one does not change the definition of the 
transverse amplitudes, see jHj for detal. 

4. One can rise the question whether it is possible to re-define the radiative decay form 
factors (for V — > 75 and other similar processes) in accordance with Eq. (j3J). For the de- 
cay V — ► -yS, we have |q| = m v — m s , so the re-definition does not affect the fact that 
F v ^s( m v, m % Q 2 — 0) depends on my and m| only [TUj . However, it is necessary to take into 
consideration that the convolution (0J) represents also a standard determination of the transi- 
tion form factor for the space-like momenta q 2 = — q 2 < 0. So, within this re-definition and 
applying q 2 = — |q| 2 = —{ray — ms) 2 one has the equality 

F v ^ lS (m 2 v , m 2 s , 0) = F v ^ lS (rn 2 v , m%, -{m v - m s f) , 

which does not look reliable. 

4 Conclusion 

We have obtained formulae, within a non-relativistic quark model, for the transitions V — ► 7S 
and V — > 7P starting from Feynman integrals. We have shown that, after a correct transition 
to the non-relativistic approximation, the decay amplitude with the emission of a real photon 
(q 2 = 0) is determined by the convolution of wave functions (see (JUT]) and where there 

is no photon-momentum dependence. This is a natural consequence of the Lorentz-invariant 
structure of the transition amplitudes, and it is a common property independent of whether 
we use relativistic or non-relativistic representations of the amplitude, see the discussion of 
formulae (3) and (4) here, as well as Eqs. (32) and (61) of paper [3]. In this point our formulae 
differ from those of the paper jB] , where there is a q-dependence in the quark wave functions. 

Another difference between papers El El and PQ 12 E] consists in the application of the 
threshold theorem j?j (see also [UJ [TJj) to the decay 0(1020) — > 7/o(980) - we discuss this 
problem in separate publication 0. 

We thank Ya.I. Azimov, L.G. Dakhno, S.S. Gershtein, Yu.S. Kalashnikova, A.K. Likhoded, 
M.A. Matveev, V.A. Markov, D.I. Melikhov, A.V. Sarantsev and W.B. von Schlippe for helpful 
and stimulating discussions of problems involved. This work is supported by the RFBR Grant 
N 0102-17861. 
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Figure 1: Quark diagram for the transition form factor V — > jS/P. 
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